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Існування та додатність розв'язків 
узагальненої моделі динаміки Гомперца 


Побудова та дослідження компартментних моделей популяційної динаміки вимагає отримання 
конструктивних умов існування та додатності розв'язків. Тоді як в лінійних випадках та таких, що до 
них зводяться, для таких задач отримані ефективні критерії, то для нелінійних систем тут виникають 
певні труднощі. У роботі розглядається узагальнена модель динаміки Гомперца, яка, на відміну від 
«скалярного» випадку, не має явного розв'язку. Застосовуючи принцип стискувальних відображень 
для такої моделі отримано конструктивні умови існування та додатності розв'язків. 


Вступ 


Динаміка Гомперца - це динаміка, яка є експоненціальною в малі моменти часу і 
прямує до деякого асимптотичного рівня у великі моменти часу | 1. 

Найпростіша модель Гомперца (запропонована Б. Гомперцом, 1779 - 1865) пред- 
ставлена рівнянням: 


де с - деяка стала. 
За формулою Гомперца - Мейкхема розраховується кількісна оцінка показників 
смертності, старіння, тривалості життя. Залежність смертності від віку |21: 


ті)- Ве". 


Модель Гомперца використовують для моделювання впливу факторів зовніш- 
нього середовища на показники здоров'я населення |31. 

Доведено, що динаміка Гомперца визначає існування фрактально-стохастичного 
дуалізму на мікроскопічному рівні надмолекулярних, клітинних систем |41. 

Моделі Гомперца знайшли практичне обгрунтування в онкології при описі росту 
ракової пухлини. Для цього до уваги беруть такі міркування. Експоненціальне зростан- 
ня - це найпростіше можливе зростання кількості клітин, що відповідає поділу клітин, 
який відбувається у рівні проміжки часу. Експоненціальне зростання не може тривати 
нескінченно. Коли розміри пухлини становлять кілька відсотків від розміру органа-мі- 
шені, орган-мішень не може повністю витримувати пухлину. В цей час експоненціальне 
зростання сповільнюється 1 прагне, аби досягти деякого асимптотичного рівня (1. 

Зв'язок між динамікою Гомперца і динамікою росту популяції клітин був вико- 
ристаний у створенні змішаної моделі, яка дає аналітичний розв'язок для клітин, що 
розвиваються, 1 тих, що перебувають в стані спокою, а також співвідношення між пара- 
метрами моделі (5). 
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У роботі (6) запропоновано узагальнену модель динаміки Гомперца з урахува- 
нням підпопуляцій пулів пухлинних клітин та їх впливу на нормальні клітини. Тому 
метою даної роботи є дослідити умови існування та додатності розв'язків моделі, 
яка може бути використана при вирішенні задач оптимального керування процесом 
лікування онкологічних захворювань. 


Основна частина 


Розглядається модель: 


У алідію -озУваій зок, 
б Зя жи З 


ге|/ оо). (1) 
з початковими умовами: 
подані, МАЛЕ М, ет; (2) 
які задовольняють: 
к ан 
Знав «МО, Ре (3) 


ізі 


До моделі (1) - (3) можуть бути зведені моделі динаміки Гомперца, що розгля- 
даються в роботі. 

У даній роботі буде встановлено умови, при яких існують додатні розв'язки (1) - 
(3), що задовольняють: 


т) з (п. (0), (є) є 5.Г1, со), | М(є) - (М. (0),..., М, (г) є Сас), (4) 


де 5Йозоо) З Повує Ск) Ук вує по 


Поряд з рівняннями (1) розглядаються їх інтегральні аналоги: 


дпа | Уалібіпчо" 


| Уп./г) 


«Урлідш-А,(НО) із їж (5) 


«дл у Хелм АМОЇ у 


Лема 1. Нехай: 
1) початкові умови системи (1) задовольняють нерівність (3); 


Кк 
2) виконується умова (До) для довільного вектора є А" :0« р «0: 
феї 


ев 
У а «У, а М. »д, а, »0, ізо. (6) 


і,ззі і,ззі 
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Тоді будь-який розв'язок (1), (2) для довільного / є (г,,со) задовольняє нерівності: 


Улпі)»о, (7) 


їзі 


м)»0, у еїп. (8) 
К 

Доведення. Припустимо, що існує перший момент часу 7, 2 1,, Коли 5 л,г,)- 0. Тоді 
ізі 


К 
внаслідок неперервності в як завгодно близькому околі г, існує /,, такий, що 5 7,(1)»0. 


їзі 


. й Мир 0. х : 
Тоді, скориставшись рівностями (5) та використавши нерівність Іп--2 В і, х0»0 іпри- 
Х 


пущення (До), отримуємо: 


б ува, «Мага 


гак із фо. ЗА (2) 


а. 


Кк кор| к 
гУліє)я У У а,п.) 
ізі 17 зі 


ай 


ох адп (ат » 0. 


і 


Отримали суперечність. Отже, (7) має місце. 
Розглянемо друге рівняння в (1). Використавши векторні позначення та формулу 
варіації сталих, маємо: 


т 


іо М) ехрід(о м Да Мо) | екріО здо М , (9) 


їо 


де 2-4)» с Зоо ік» ехр( ) - матрична експонента, Іл М(г) - (іп М, (с). а М, (г))". 


Звідси бачимо, що з умов (3) випливає (3). 
Теорема 1. Нехай виконано умови леми 1. Тоді існує єдиний розв'язок задачі 


0)- 0). 


Доведення. Зауважимо, що розв'язок рівнянь для М, (г) існує і є єдиним, як 
тільки такі властивості має розв'язок рівняння для /, (г). Для цього слід скористатися 


формулою (9). 
В силу леми 1 розв'язок системи (1) повинен задовольняти нерівності (7), (3). 


Тому, використавши нерівність Іп о і, х,0»0 в (1) отримаємо: 
Х Х 


аа з ал м )«ак ел лек (10) 
й Ух 


з-Ї 
Отже, система (1) - (3) мажорується лінійною системою з тими ж початковими умо- 
вами, розв'язок якої неперервний 1 існує при всіх / 2 І,. Неперервність правих час- 


тин (1) 1 їх похідних в околі початкових умов гарантує локальне існування і єдиність 
розв'язку |7). 
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Наслідок 1. Нехай виконуються умови теореми І та виконується умова (А!) - 
або матриця 8 -|,| вироджена, або матриця В - не вироджена 1 для довільного 


ізаі К 


К 
вектора 5 - (2 б) є Кк 323 2 д виконується нерівність: 


ізі 


К 


р лови 


о Уа, адбдачьк (п 


їз 
З. Є іззеі 
Р 


Утпі)«в. (2) 


Доведення. У випадку, якщо матриця В - вироджена, то розв'язки, для яких 


Кк 
Улп()є8, є точками рівноваги і нерівність (12) виконується. 


їзі 


Розглянемо випадок, коли В - невироджена. Тоді з рівняння (1) видно, що 


К 
жодна точка, для якої У 7, (г)- 9 , не може бути точкою рівноваги. 


ізі 
4 
Припустимо противне, тобто нехай існує момент часу 1, коли 5 п(/)-0. 
ізі 


Оскільки точка (г (г)... (г) не є точкою рівноваги, то в як завгодно близькому околі 


Кк 
існує момент часу 7,, коли Упі)» 9. Тоді, використавши рівняння (1) і нерівність 
ізі 


то-«о--і х,2»0 маємо: 
Х 


Отримали суперечність. 
Наслідок 2. При виконанні умов теореми 1 оператори А. (п), ізі,К та 


А, (те), М(г), )-іп, що означені на відповідних банахових просторах, задо- 


вольняють принципу стискаючих відображень, а їх нерухомі точки є розв'язками 
задачі (1) - (3), які задовольняють умови (7), (8). 
Висновок 


Отже, в роботі розглянуто модель узагальненої динаміки Гомперца. Дана модель 
враховує різні пули пухлинних клітин та підпопуляції нормальних клітин. Запропоно- 
вано конструктивні умови існування, єдиності та додатності розв'язків моделі. 
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Существование и положительность решений обобщенной модели динамики Гомперца 
Построение и исследование компартментньїх моделей популяционной динамики требуєт получения 
конструктивньх условий существования и положительности решений. Тогда как в линейньх случаях 
и таких, которью к ним сводятся, для таких задач полученьт зффективнье критериий, то для нелинейньтх 
систем здесь возникают определеннье трудности. В работе рассматриваєтся обобщенная модель 
динамики Гомперца, которая, в отличие от «скалярного» случая, не имеет явного решения. Применяя 
принцип сжимающих отображений для такой модели полученьт конструктивнь6е условия существования и 
положительности решений. 


У.Р. Магізепіуик, І.Е. АпагизИиснаєї, І.5. СуосуеїзКа 

Ехі5(епсе апа Розійуїгу ої 501ийоп5 ої Сепегаїгед Модсеі ої Сотрегі7/5 рупатіс5 

Виійпе апа іпуезісайоп ої сотрагітепі подеїя ої рорціайоп Фупатіся гедшіге соп5ігисіїує сопаїіопе 
Гог ехізіепсе апа розійутйу ої зоїийопя. М/Пе іп Ппеаг са5е5 апа іп са5е5 гаїзед 10 пет ейесйує сгіегіа 
аге Гоцпа, Феп Феге аге 5опле ФіНбсиіце5 їп попіїпеаг зузієтя. ТРе агіїсіє деа5 мії Фе репегайоей подає 
ої Сотрегі7'8 дДупатіся, уубісі 15 діНегепі гоп "5саіаг" сазсапа доезп'ї рауе ехрісії 5оїийоп. Ву цз5іпе 
ргіпсіріе ої сопігасіїпо пзаарріпе5 Гог Бі5 плоді Фе соп5ігасбує сопайопз бог ехізіепсе апа розійутсу ої 
50їийоп5 Баме Бееп Гоцпа. 


Стаття надійшла до редакції 15.03.2011. 


«Штучний інтелект» 272011 163 


